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ABSTRACT 



The treatise seeks to find the optimum shape of a thin 
beam under the influence of lateral buckling. The specific 
formulation is by means of the isoperimetric problem of the 
calculus of variations. An energy method approach yields 
the governing equations for the problem. An analytic 
solution was not obtained due to the nonlinearity and 
coupling conditions of the equations. 

A variable height lateral buckling problem is formulated 
as an alternate attempt to find the optimum design. Through 
Bessel equation analysis numerous designs are obtained and 
the resulting buckling loads are calculated. The largest 
buckling load corresponds to a beam design which appears to 
be very close to the optimum shape. 



2 



TABLE OP CONTENTS 



I. INTRODUCTION 9 

A. THE OPTIMIZATION PROBLEM 9 

B. SCOPE OF THE THESIS 13 

II. PROBLEMS OP OPTIMAL STRUCTURAL DESIGN 13 

A. THE BASIC PROBLEM OP STRUCTURAL MECHANICS 13 

B. CALCULUS OP VARIATIONS 14 

C. THE ISOPERIMETRIC PROBLEM 15 

D. POTENTIAL ENERGY 16 

E. PRINCIPLE OP MINIMUM POTENTIAL 17 

F. A STRUCTURAL OPTIMIZATION PROBLEM 18 

III. . THE LATERAL BUCKLING PROBLEM 22 

A. GENERAL DESCRIPTION 22 

B. DEVELOPMENT OF CLASSICAL SOLUTION 23 

IV. A VARIABLE HEIGHT BUCKLING PROBLEM 32 

A. DESCRIPTION OF THE PROBLEM 32 

1. Non-integer Value v 35 

2. Integer Value of v 37 

B. CALCULATION OF AMMENDED LOADS 38 

C. EXTENSION OP THE PROBLEM 43 

D. LATERAL BUCKLING OP A CANTELEVER 45 

E. COMPARATIVE ANALYSIS OF SOLUTION 47 

F. VALUES OP n GREATER THAN ONE 52 

V. CONCLUSIONS 54 

APPENDIX A: A STRUCTURAL OPTIMIZATION PROBLEM 56 



3 



APPENDIX B: EFFECT OF THE LOAD P AT A DISTANCE 

FROM THE CENTROID 62 

APPENDIX C: TRANSFORMATION OF EQUATION (4.11) 

TO A BESSEL EQUATION 64 

APPENDIX D: BESSEL EQUATION CURVES 66 

APPENDIX E: VARIABLE HEIGHT DESIGNS 76 

COMPUTER PROGRAM 108 

LIST OF REFERENCES 12 3 

INITIAL DISTRIBUTION LIST 125 

FORM DD 1473 126 



4 



LIST OP FIGURES 



1. Column Shapes Analyzed 18 

2. Buckled Beam 22 

3. Beam Loaded at the Centroid 23 

4. Lateral Buckling Configuration 25 

5. Variable Height Beam 29 

6. Variable Height Beam 32 

7. Convergence of Power Law 39 

8. Variable Height Beam (n=0.l4) 42 

9. Graph P* vs n (S.S. beam - case 1) 4 3 

10. Variable Height Eeam 44 

11. Graph P* vs n (S.S. beam — case 2) 46 

12. Graph P* vs n (cantilever beam ) 48 

13* Variable Height Cantilever 49 

14. Moment Configuration 51 

15. Variable Width Beam 56 

16. Optimum Shape Variable Width Design 6l 

17. Effect of P at the Top of the Beam 63 

18. -27. Bessel Equation Curves 66-75 

28.-59. Variable Height Eeam Designs 76-107 



5 



LIST OF SYMBOLS 



A 



a 

b 



b 

o 



D 



E 

F 



f 

G 



h 

h 

< 

I 

i 

J 

j 

K 



k 

L ,1 

M 

N 



n 

P 



cross sectional area 
constant of integration 
variable width 
constant width 
material constant 
material constant 
domain 

Young's modulus 

functional 

function 

shear modulus 

variable height 

constant height 

integral function 

vector subscript 

polar moment of inertia 

vector subscript 

dimensionless constant 

dependent variable subscript 

length 

moment 

column load 

exponent for design 

concentrated load 

uniform load 



6 



s 

T 

U 

V 

V 

o 

U,VjW 

x 3 y,z 

a 

3 

V 
6 
e 

X 

v 

P 

a 

4 > 



surface 

total potential energy 
strain energy 

potential energy external forces 
volume 

displacements 

coordinate axes 

constant of Integration 

Bessel equation parameter 

(1 - 2x/l) 

variation 

strain 

principle axes 
Lagrangian multiplier - 
Bessel function order 
Bessel equation argument 
stress 

specific strain energy 



7 



ACKNOWLEDGEMENTS 



The author wishes to express his appreciation to his 
advisor Associate Professor David Salinas for his construc- 
tive supervision and enthusiastic encouragement throughout 
this work. He also wishes to thank Associate Professor 
Craig Comstock for his help in overcoming mathematical 
difficulties encountered during the thesis. 

The Naval Postgraduate School Computer Center provided 
facilities for computer work. 



8 



I. INTRODUCTION 



A. THE OPTIMIZATION PROBLEM 

The performance of an elastic body is measured by the 
maximum load which the body is able to withstand. The 
optimization problem considered in this treatise is that 
of determining the specific design of an elastic structure 
such that the structure will exhibit a maximum performance 
capacity. Herein, the term design refers to the actual 
physical distribution of material within an elastic body. 

The optimization problem is formulated by the calculus 
of variations. This method offers certain advantages in 
the development of a theory for optimization of an elastic 
body. The specific formulation of the problem is by means 
of the isoperimetric problem of the calculus of variations. 
By this method the problem becomes an extremization of the 
total potential energy of the system with a superimposed 
volume constraint. The extremization of the functional, 
augmented by Lagrangian multipliers, leads to a set of 
Euler equations which govern the design of the structure. 

The calculus of variations is no new technique for 
obtaining an optimum design. Lagrange [1] formulated an 
optimum column problem by this means in 1773* He sought to 
determine the shape of a column for which the buckling load 
was largest among all columns of given length and volume. 
Although he failed in obtaining the correct solution, his 
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recognition of the problem belonging to the calculus of 
variations inspired others to continue research in this 
field. Clausen [2] later provided a solution to this prob- 
lem for columns of circular cross sections. Keller [3] 
recently showed that the strongest column is not one of a 
circular cross section, but instead, is one which has an 
equilateral triangle as a cross section. 

The isoperimetric problem for optimum design received 
an additional stimulus from Keller and Tadjbakbsh [ 4 ] . This 
work provided a variational technique for minimum volume, 
elastic design of columns. Taylor [5] treated the same 
problem by means of an energy method approach. In his work 
Taylor showed that the same governing equations may be 
obtained by direct extremizati-on of a functional, which 
consists of the system potential energy and a volume con- 
straint. Taylor's approach is the essential method used 
for optimization in this thesis. The method is particularly 
useful since it provides a simple and direct means for iso- 
perimetric problem formulation. 

An extension of the Taylor approach was accomplished by 
Salinas [6] and provides some conditions for the validity 
of the potential energy formulation. Additionally, the 
author demonstrates the equivalence of maximum load, minimum 
volume optimum designs and presents a general formulation for 
optimization of elastic structures. 
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B. SCOPE OF THE THESIS 



Using the Taylor approach, the optimization problem is 
discussed and developed by use of Lagrangian multipliers 
and the calculus of variations. The resulting set of Euler 
equations are seen to be the system equilibrium equations 
and an optimality condition. Example problems are included 
to illustrate the usefulness of the Taylor method. These 
include a strength problem for a uniformly loaded beam and 
an optimum design problem for a particular case of column 
buckling . 

The method is then extended to the case of a thin beam 
subjected to lateral buckling. Under lateral buckling, the 
beam is considered to buckle in a plane perpendicular to 
the vertical plane of loading.- The buckling configuration 
consists of rotation and sidewise bending away from the 
horizontal axis of symmetry. The resulting equations 
governing an optimum design are found to be non-linear. In 
addition, the equations are coupled in rotation and bending. 
All attempts for an analytic solution proved futile. 

An attempt to find the optimum shape has been made by 
assuming geometrical conditions for instability and design 
of the beam. A variable height problem has been formulated 
by means of a method described by Federhofer [71 . Numerous 
designs are obtained by use of Bessel equation analysis and 
a predicted optimum shape is forseen. 

Throughout the discussion all problems are considered 
only for the elastic range of the material. Additionally, 
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the material is considered to be homogenous and isotropic. 
All problems analyzed are restricted to rectangular cross 
sections . 
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II. PROBLEMS OF OPTIMAL STRUCTURAL DESIGN 



A. THE BASIC PROBLEM OF STRUCTURAL MECHANICS 

The classical theory of structural mechanics is concerned 

with the behavior of continuous elastic bodies which are 
subjected to specific loading conditions. In general the 
components which define mechanical behavior are the displace- 
ment vector u, and the stress and strain tensors, a., and 
* iJ 

£ ij* 

The strain tensor is explicitly defined by the geometry 
of deformation of an elastic body. Hence, the strain- 
displacement relation may be represented by 

£ = £ (u) - (2.1) 

From the constitutive laws of elasticity it is known that 
stress is related to strain, that is 

a = a(e) (2.2) 

Thus, stress and strain are both representable as functions 
of the displacement. In classical theory any one of the 
three quantities, u, a, and e is sufficient for determining 
the other two. 

The mechanical state of a structure is defined by a 
boundary value problem which is comprised of equilibrium 
equations and the associated boundary conditions. Thus, a 
classical problem in structural mechanics consists of using 
the constitutive relations between stress, strain and 
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displacement to derive and solve the governing equilibrium 
equations and their associated boundary conditions. 

B. CALCULUS OF VARIATIONS 

Problems concerned with the determination of extreme 
values of integrals whose integrands contain unknown func- 
tions below to that field known as the ’’calculus of 
variations" [8] . Such an extremal refers to the maximum, 
minimum or stationary value of a functional. A functional 
is a real valued function of a function. 

Only the integral functional 

I[f(x)] = / F[f(x)]dx (2.3) 

D 

will be considered herein, where D is the n-dimensional 
domain over which the n-dimensional vector x ranges. The 
function which extremizes I[f(x)] is denoted as f#(x) and 
the associated value of the functional is I* (i.e. I* = 
I[f*(x)]) . 

In general, problems in structural mechanics concerning 
optimization will assume the form 

/ F(x;u,u' ,u" ,4> )dx (2.4) 

R 1 

for a one dimensional case. A one dimensional case is 
generally characterized as having one dimension, its length, 
much larger than the other dimensions. In the above form x 
represents an independent variable. The quantities u, 4> and 
their respective derivatives represent dependent variables. 
In this thesis u and 0 represent displacement functions or 



state variables. The nature of the specific problem governs 
the number of dependent and independent variables. 

Performing a variation with respect to a functional is 
denoted by the symbol 6. From [8] it can be shown that a 
stationary value of the functional occurs for 6F = 0 . This 
operation yields what is known as a set of Euler equations. 
For example, if equation (2.4) is represented by the par- 
ticular form 



I = / F(x ;u,u' ,u" ,<f> ,<j> ' ,A)dx (2.5) 

D 

then performing the extremal operation yields three Euler 
equations [91 • 
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( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 



plus a set of natural boundary conditions. 



C. THE I SOPERIMETRI C PROBLEM 

The isoperimetric problem in the calculus of variations 
provides for extremizing functionals whose dependent vari- 
ables are constrained by an integral relation. The integral 
constraint is accommodated by constant Lagrangian 
multipliers . 

As an example, suppose that the maximum (or minimum) 
value of an integral 
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\ 



02.9) 



I = / F(x ;u ,u ’ )dx 
R 

is sought, subject to the condition that another integral 

J = / G(x;u,u’)dx (2.10) 

R 

is known to have a constant value. Using the method of 
Lagrange multipliers [8] an augmented functional 

I+AJ = / [F ( x ;U ,u ' ) + G(x;u,u')]dx (2.11) 

R 

may be constructed and its free extremum considered. The 
necessary condition for a stationary value is the Euler 
equation 

3 (Ft 0) _ 3 (F+ G) n (0 

3u 3x 3u» Kd.id) 

where, F has been replaced by F-* = (F+AG). 

D. POTENTIAL ENERGY 

The total potential energy of a conservative elastic 
system can be expressed as 

T = U - V . (2.13) 

U represents the total strain energy of the structure, and 
V is the potential energy of external forces associated with 
the displacement of the system under any set of given load- 
ing conditions. The total strain energy of the system is 
the volume integral 

U = \ f a. .e. .dV (2.14) 

2 y 1J 

where a., and e are components of the stress and strain 
X J ij 

tensors . 
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The potential energy of external forces is defined as 
the inner product of the external force vector P and the 
associated surface displacement vector u, that is, 

V = / P * u ds 
s 

where s is the domain over which P is defined. 

E. PRINCIPLE OF MINIMUM POTENTIAL ENERGY 

A powerful tool in the field of structural mechanics is 
the principle of minimum potential energy. This theorem 
states that the total potential energy of a conservative 
elastic system attains a minimum, with respect to variation 
of kinemati-cally admissible displacement functions, when 
the displacement function satisfies equilibrium. 

The total potential energy has previously been defined 
in terms of strain energy and external work. These quanti- 
ties may be expressed as functions of the displacement 
vector u. Thus, the equilibrium equations and associated 
boundary conditions are derivable from a variation with 
respect to the displacements of a potential energy functional . 
The variational operation 

<5 U (U-V) = 0 (2.16) 

can then be performed for which the resulting Euler equations 
physically represent the system equilibrium equations. The 
optimization method considered in the thesis is based on 
this important principle. 
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F. A STRUCTURAL OPTIMIZATION PROBLEM 



As an example of structural optimization, consider the 
problem of obtaining the optimum shape for column buckling. 
The problem may be characterized by the statement : "For a 

column of given length and volume of material, determine the 
column shape for which the Euler buckling load is a maximum. 

The column shape is considered to be such that all cross 
sections are similar. The moment of inertia I(x) is 
assumed to be related to area, A(x), that is 

El ( x ) = kA n (x) (2.17) 



where E represents Young's modulus and k is a constant 
depending on the section considered. For example, the 
moments of inertia of rectangular sections of variable- 
width-constant height, and variable height-constant width 
are (Fig. 1), 




Figure 1. Column Shapes Analyzed, (a) variable width - 
constant height, (b) variable height - 
constant width. 
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(n = 1) 



( 2 . 18 ) 



El 



yy 



Eh 3 b(x) 

12 



kA ( x ) 



for case (a) and 



EI yy = " bh i 2 X) = k A 3(x) (n = 3) (2.19) 

for case (b ) . 

For the Euler buckling problem the total potential 
energy Is 

£ 

T = \ / (EIv" 2 - Nv ' 2 ) dx (2.20) 

^ 0 

where v represents lateral displacement for the buckled 
state. Introducing an isoperimetric condition for a con- 
stant volume of material V gives 

o 

£ 

/ A(x)dx = V (2.21) 

/-v o 



The volume constraint is accommodated through the 
formation of the functional T* , defined by 

T* = T - XV (2.22) 

o 

where X is a constant Lagrangian multiplier. Combining 
(2.17), (2.20) and (2.22) the constrained functional is 

l 

T* = / (- A n v" 2 - | v' 2 - AA)dx . (2.23) 

0 2 2 

The governing equations, in addition to the isoperimetric 
equation, are obtained by performing a variation with 
respect to v and A to obtain 
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0 



(2.24) 



6 y T* = 0: (kA n v" ) " + Nv" = 

6 a T* = 0: |nkA n_1 v" 2 - A = 0 (2.25) 

Equation (2.24) is the equilibrium equation and (2.25) an 
optimality condition. 

The non-linear equations (2.24), (2.25) and (2.21) define 
the optimum design with respect to constant volume - maximum 
load. In general these are extremely difficult to solve 
in closed form. However, an exact solution can be obtained 
for case (a) where n = 1. These governing equations then 
become 

(kAv")" - Nv" = 0 (2.26) 

\ kv" 2 = X (2.27) 

i 

f A(x)dx = V q (2.21) 



Equation (2.27) implies constant 
v" = /2A/k 

Substituting (2.28) into (2.26) 

A = ic ^ ~2 + C 1 X + C 2 ^ 



bending curvature, i 



and integrating twice 



e . 



( 2 . 28 ) 

gives 



(2.29) 



Similarly 



v = /2A/k + C 3 x + C^) 



( 2 . 30 ) 



If the column is considered to be simply supported, the 
boundary conditions for displacement are 



v(0) = v(£) = 0 



(2.31) 
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for which c ij = 0, and 



v ( x ) = /2X/k 



The moment conditions are 

M(o) = M(£) = 0 
£ 

from which = 0, and 

.2 



A (x ) = | (2L. 



£x 



) . 



( 2 . 32 ) 



(2.33) 



(2.3*0 



The critical load is now determined from the isoperimetric 
condition 



£ 

/ 

0 



N ,x£ _ 

k 1 2 2 ' 



V 

o 



(2.35) 



for which 



N 



_12kV 0 -Eh V 



o 



optimum 



( 2 . 36 ) 



Substituting this expression into (2.3*0 determines the 
optimum design in terms of known constants, that is 



A(x) = — (£x - x 2 ) . 
12 £ ^ 



(2.37) 



Note that the Lagrange multiplier does not appear in either 
of the expressions for A(x) or N. This is in distinction 
to strength optimization problems where X is a representa- 
tion of maximum strength and remains in the final analysis 
(see Appendix A). 



21 



III. THE LATERAL BUCKLING PROBLEM 



A. GENERAL DESCRIPTION 

The problem of lateral buckling is important in the 
design of beams which have no lateral support and possess a 
shape such that the beam height is large in comparison to 
the width. Such a configuration may become laterally 
unstable when the loading reaches a critical value. The 
resulting failure or "buckle" is a combination of lateral 
bending and twisting away from the horizontal axis (Fig. 2). 
Prior to buckling the load needed to cause this failure may 
be considerably less than the corresponding load that would 




Figure 2. Simply Supported Beam in a Lateral Buckled State. 
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cause bending failure. This is due to a large flexural 
rigidity of the beam in the vertical plane of bending as 
compared to a much smaller lateral bending resistance in 
the plane perpendicular to the applied load. 

A brief but informative discussion of the historical 
development of the lateral buckling problem is given by 
Bleich [12] and will not be included here. The methods 
used for development of the problem in this thesis are sim- 
ilar to those used by Sechler [13] and Timoshenko [1*1] . 

B. DEVELOPMENT OF CLASSICAL SOLUTION 

While the solution to the classical lateral buckling 
problem may be found in a number of texts, its inclusion 
herein is considered necessary for understanding the 
descriptive geometry of the problem and later development 
of the problem of a variable height solution. The example 
to be considered is that of a simply supported beam loaded 
at its centroid by a concentrated load P at the origin 
(Fig. 3). The centroidal loading condition is necessary to 
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simplify the problem. Loading at any point other than the 
centroid introduces an additional work term which leads to 
some mathematical difficulties. This condition is described 
in Appendix B. 

It is assumed that the load P remains in the vertical 
xy plane throughout buckling. The ends of the beam are so 
constrained as to prevent rotation in the yz plane. The 
buckled configuration is shown in Figure 4. 

The deflection of the beam is defined at any point x 
from the origin by the components u, v and <J> . Components 
u and v represent the displacement of the centroid in the 
z and y directions while <J> is the angle of rotation of the 
cross section in the yz plane. Thus, the boundary condi- 
tions, for the problem are 

uU/2) = <J>U/2) = 0 (3.1) 

- u ' ( 0) = <f>'(0) = 0 (3.2) 

The fixed coordinate axes are labeled x, y and z. In 
the buckled state the principle axes of symmetry at any 
section mn are defined as £, n, and £ . The usual right 
hand rule convention for vectors is used to describe posi- 
tive moments with respect to the described axes. 

By the usual assumptions of small deflection theory the 
equations for principle moments can be written 
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( 3 . 3b ) 



H 2 r 

El ^-§- a El u" = -M 
n dC 2 n n 



GJ || a GJ<t> T = 



(3-3c) 



These moments are physically described as follows: 



i) M is the usual bending moment due to the 
5 load P 

ii) M is the bending moment associated with 
n lateral buckling 

iii) M- is the twisting moment associated with 
^ lateral buckling. 

With respect to the fixed x, y and z axes the above 
moments can be represented [131 as 



m 5 


= M + 

X 


M u* 
z 


= GJ<j> * 


(3.4a) 


M 


= M + 


M <J> 


= -El u" 


(3.4b) 


n 


y 


z 


n 


M 


= -M <J> 


+ M 


= EI_v" 


(3.4c) 


C 


y 


z 


C 


Let 


= GJ and C 2 


= El . Differentiating, 
n 


the first 



equation with respect to x and substituting for u" from the 
second results in a general equation which can be written 



M 



<p" + 



C 1 C 2 



, , MM 

<t> = — m* + — m'u* - — y 

9 C 1 x C 1 z u C 1 C 2 



(3.5) 



For the problem considered let the deflection at the 
beam center (maximum deflection) in the z direction be 
denoted by u^. The moments at any distance x are determined 
as 



M x = -P/2 (u 1 -u) (3.6a) 

M = 0 (3.6b) 

y 
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(3.6c) 



M = P/2 U/2-x) 

Zi 



Differentiating and substituting into equation (3*5) 
yields 



4>" 



^2 ( i 
P (o “ 







0 . 



(3.7) 



Introducing the new variable 6 = (i/2-x) and letting 
2 

2 P 

k = {.p p equation (3.7) becomes 
40 1°2 

4 >" + k 2 B 2 <t> = 0 . ( 3 . 8 ) 

Equation (3-8) is now in the form of a Bessel equation for 
which the general solution has the form 

4> = B h ^^(l 8 2 ) + a 2 J_ 2< (| 6 2 )] (3-9) 

where J, and J 2 are Bessel functions of the first kind of 

4 

orders h and respectively [ 1 4 ] . For the simply supported 
beam the boundary conditions with respect to <f> are 



<t> = 0 at 


6 = 


0 


(i.e. x = i/2) 


(3- 


10a) 


<t> ' = 0 at 


6 = 


i/2 


(i . e . x = 0) . 


(3- 


10b) 



The first boundary condition yields a 2 = 0 . Differentiation 
of (p with respect to 6 now yields 



d£ 

d6 



a x kB 



3/2 



J. 3/1) ( | e 2 ) 



(3.11) 



Applying the second boundary condition shows that buckling 
occurs when 



•J 



-3/4 



(- T _) = 0 



( 3 . 12 ) 
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From a table of Bessel functions [15] the first root of 
(3.12) is 



kl 



2 



8 



1.058508 . 



(3.13) 



Substituting for K, the buckling load is found to be 



P 

cr 



16.936 



(3.1*0 



Equation (3*1*0 represents the exact solution for the prob- 
lem considered and will be used as a comparative value for 
critical loads calculated for variable cross sectional 
designs later. 



C. OPTIMIZATION TECHNIQUE APPLIED TO LATERAL BUCKLING 



The energy formulation described in Section II is now 
used to develop the governing equations for an optimum shape 
of a thin beam which undergoes lateral buckling. The 
assumptions made during the preceeding section apply. Addi- 
tional assumptions include neglecting the strain energy due 
to "warping" [14] and elongation of the beam. 

The case considered is that of a variable height-constant 
width beam (Fig. 5). The principle moments of inertia with 
respect to the buckled state are 



I 

n 



= T? b o h(x) 



O 

12 



A ( x ) 



J ' 3 



b" 5 h(x) = 



o 



o 

3 



A(x) 



(3.15) 

(3.16) 



The strain energy due to lateral bending may be written 
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U 1 = ; 
1 0 



1/2 



C 1 Au" 2 dx 



(3-17) 



where C, = Eb /12 . The strain energy associated with 
1 o 

twisting of the structure is 



1/2 



U ? = / CpA<t>’ 2 dx 
0 



(3.18) 



where = Gb^/3. 

As the beam buckles the load P is lowered an amount which 
may be represented as [14] 

<t>u" U/2 - x) dx . (3.19) 

The potential energy of external forces associated with the 
buckling configuration is determined by 



1/2 



V = P / 
0 



<J>u" U/2 - x) dx . 



( 3 . 20 ) 



The total potential energy of the system is U - V or 



1/2 



T = / [C, Au" 2 + CpA<j> ' c - P^u" ( 1/2 - x) ] dx. (3-21) 

0 1 ^ 



,2 



To formulate a functional by the method described earlier. 
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an isoperimetric constraint is introduced for a given volume 



of material. 



1/2 

f Adx 
0 



= V 



( 3 . 22 ) 



The augmented functional may be written 

A/2 P P o 

T* = / [C-. Au" + C A<j>' - P4>u"(^-x) - AA] dx. (3.23) 

0 1 

Performing a variation with respect to the variables u, <j> 
and A yields three Euler equations : 



6 T# = 0: 
u 


[2C 1 Au" - P4> (j 


- x)] 


= 0 


(3.24) 


6.T* = 0: 
<P 


2 C 2 ( A4> ' ) ' + Pu" 


(* - 
' 2 


x) = 0 


(3.25) 


O? 

> 

i-3 

He 

II 

O 


C 1 U " 2 + C 2 <t> ' 2 ~ 


A = 


0 . 


(3.26) 



Note that the variation with respect to the state variables 
u and yields equilibrium equations. Variation with 
respect to the design variable A introduces an optimality 
equation. The associated boundary conditions are: 

u(f) = <*>(§) = 0 (3.27) 

u* (0) = <p ' (0) = 0 (3.28) 

Integration of equation (3.24) twice and evaluating at the 
boundary conditions determines 

2C 1 Au" - P<f> ( | - x) = 0 (3-29) 

Integration of (3*25) yields the result 

C 2 A<4> t + u* (j - x) + (u-j^-u) = 0 (3.30) 
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where represents the maximum deflection at the beam 
midspan . 

The above equilibrium and optimality equations along 
with the isoperimetric condition (3*22) define the optimum 
shape for the structure. These governing equations are noted 
to be non-linear and an attempt to solve them has been 
unsuccessful. Unfortunately, these equations resist uncoupl- 
ing, which precludes getting to the isoperimetric condition 
as was accomplished easily in Section III. Additionally, no 
computer techniques were found which could solve systems of 
simultaneous nonlinear equations. 

After considerable efforts to solve the above non-linear 
system proved fruitless, an attempt to formulate a variable 
height solution from a specified polynomial was undertaken 
in hopes of predicting or converging to an optimum solution. 
The method used and results obtained are included in the 
following section. 
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IV. A VARIABLE HEIGHT BUCKLING PROBLEM 



A. DESCRIPTION OP THE PROBLEM 

The problem of lateral buckling of a variable height, 
constant width beam is now considered. Again, the problem 
to be investigated is that of a simply supported beam, 
loaded at the midspan centroid. However, the height now 
varies according to the power law 



where h Q is the height at the center section of the beam 
(Fig. 6). The value of n is restricted to lie between the 
interval 0 £ n < 1 for reasons which become physically 
obvious in later analysis of the problem. The idea for 
using this type- of formulation was based on the results 




(4.1) 



y 





X 



Figure 6. Variable Height Beam Based on the 

Power Law h (l-2x/Jt) n . 
o 
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obtained for a variable height cantilever problem by 
Federhofer [7] . 

Differentiating equation (3* 4a) with respect to x and 
noting that the previous constants and C 2 are now vari- 
able functions of x, reduces (3.4a) to the form 

C,<|>" + C _*<!>» = M* + M’u r + M u" (4.2) 

1 1 x z z 

Substituting (3. 4b) into (4.2) for u” and recalling the 
moment conditions from (3-6) (and their respective deriv- 
atives) reduces (4.2) to 



P 2 (- - x) 2 

4>” + ^ C'<j>’ + —fSr-n 4> 






= 0 



(4.3) 



Letting 



/ , 2x, 

Y = a - T ) 



(4.4) 



and using the relations 
C 1 = Eb 3 h(x) 

^2 = \ Gb ^ 

the lateral and flexural stiffness at the beam center may 
be designated as C° and C°. Thus, the new variable stiff- 
nesses can be represented by 



c i = c l (i - ¥> n = c ^ n 



C 2 - C° <1 - f) n = C°Y n 



(4.5a) 

(4.5b) 



Using the chain rule 



33 



d<j> _ clx 

dy ~ dx dy 



(4.6) 



and noting 



x = | (1-y) 



A .. = O. = O. 
v 2 2 2 

dx l dy 



(4.7) 



Similarly 



dC. 



JL 111 d^_ _ 4_ n d^_ 

dx dx ^2 y dy 



(4.8) 



Substitution of (4.7) and (4.8) into equation (4.3) yields 



,2. .. p 2 0 4 2(l-n) 

d 4> + n d^_ + P % y 

dy 2 Y dY 64C°C° 



= 0 



(4.9) 



The dimensionless parameter K may be defined 



PV = k 2 



6 4C°C° 



(4.10) 



Equation (4.9) may now be written 



+ K 2 Y 2(1 - n) t - 0 . 

dy ' Y 



(4.11) 



Transformation of equation (4.11) into a Bessel equation can 
be accomplished by the following substitutions. 



P = 



K 2-n 



Y 



9 = y 



2-n 
1-n 

2 B (y) 



(4.12a) 

(4.12b) 



Equation (4.11) is then reduced to the form 
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(4.13) 



*4 ♦ i - * a - 4> 

dp p dp p 



= 0 



where v is representative of 



\j - I"* 1 . 

V " 2 ( 2-n ) ' 



(4.13a) 



Verification of the transformation is accomplished in 
Appendix C. 

1 . Non-integer Value of v 

For the non-integer values of the parameter v 
between zero and one equation (4.13) has the following solu- 
tion [16] 

6 = a 1 J v (p) + « 2 J_ V (P) (4.14) 



where J v and J_ v are Bessel functions of the first kind and 
order v. 

1-n 

Recalling (4.12b), (i.e. <p = Y ^ 8(y)), the complete 

solution is determined by 



1-n 

4> = Y 2 [a 1 J v (p) + a 2 J_ v (p)] 



(4.15) 



The boundary conditions associated with the lateral buckling 
problem are: 

<f >( y = 0) = 0 (i.e. x = —) (4.16a) 

(Y = 1) = 0 (i.e. x = 0) . (4.16b) 

From the first condition and expansion of J_ v , it is con- 
cluded that a 2 = 0 . Therefore (4.15) becomes 
1-n 

<J> = Y 2 [a 1 J v ( p)] (4.17) 
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for which differentiation yields 



d<J> 

— = a, 
dy 1 



- 1 -n 1 -n 

1 “ n 2 t ^ ^ i 2 Tt/\ dp 

— y J v (p) + y J^(p) 57 



It is additionally known that 



(4.18) 



dp 

dy 



^7= Ky 



1 -n 



(4.19) 



and [ 17 ] 



j» (p ) = _ j (p) + j , 

v K v +1 K p v(p) 



(4.20) 



Evaluation of at the second boundary condition with 
y = 1 and p = K/(2-n) yields the solution 



J.. ^ J. 



v v 2 -n 1 -n v +1 2 -n 



(A) 



(4.21) 



The buckling loads may now be determined by solution of 
equation (4.21). This may be accomplished by independently 
calculating each case 



J v - 0 



(4.21a) 



and 



K 



(A) - 0 



(4.21b) 



1 -n v +1 2 -n 

and determining the initial intersection of the resulting 
Bessel equation curves. For each considered value of n, K 
is found numerically from the argument of J. A computer 



* 

Handbook of Mathematical Functions, National Bureau of 
Standards Applied Mathematics, Series 55 , P* 363 , June 
1964 . 
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solution for this problem is included in the Computer Pro- 
gram section of the thesis. The Bessel equation curves 
which show the intersecting roots are included in Appendix D. 

When n=0 the beam design reduces to that of the rectang- 
ular case. The computed solution from equation (4.21a) and 
(4.21b) provides an argument value of I.O 585 I. Thus, the 
critical load is determined by 

7T- = 1.05851 (4.22) 



or 

16.936 Jc^c° 2 

P cr = 



(4.23) 



precisely the value determined previously. The values for 
other solutions are discussed in a later section. 

2 . Integer Value of v 

When the value of n is equal to one the resulting 
value of v becomes zero and equation (4.13) is reduced to 
the form 

+ i ^+B = 0. (4.211) 

dp 2 p dc 



The complete solution of this equation can be shown to be 
of the form [17] 

6 = a 1 J Q (p) + a 2 y o (p) (4.25) 



where y is known as a Weber function or Bessel function of 
' o 

the second kind. Substituting for 6 , the above equation is 
transformed to 
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1-n 

<P = Y 2 ta 1 J Q (p) + a 2 y Q (p)] . (4.26) 

This solution is valid only if v is of an integer order. 
Thus, in the restricted range n is allowed only to be equal 
to one, and the y multiplier must be to the zero power. 
Equation (4.24) reduces to 

4> = a 1 J Q (Ky) + a 2 y o (Ky) . (4.27) 

0 

The value of y Q (Ky) for y = 0 (i.e. x = is which 
implies that is equal to zero. Differentiating <t> with 
respect to the remaining term and evaluating at the second 
boundary condition (y = 1, x = 0) requires that 

-J X (K) = 0 (4.28) 

for the solution to (4.26). This value is found C 1 81 to be 
either the trivial case (i.e. zero) or 3-8317- The latter 
result implies a singularity condition with respect to the 
previous non-integer solutions. A graphical representation 
(Fig. 7) of values of n versus the roots of equation (4.21) 
does, in fact, imply the resulting value converges to zero 
as n approaches one as a limit . Therefore , the trivial case 
is the solution. The beam is apparently unstable for any 
load when in this configuration. 

B. CALCULATION OF AMENDED LOADS 

Extension of the problem is now made in an attempt to 
find the design among all those analyzed that maximizes the 
buckling load. Initially h Q was described as a fixed height 
at the beam midspan. This center height is now allowed to 
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X 




VALUE OF N 

Figure 7- Convergence of the Root x as N approaches One. 



vary such that each variable design studied will have the 
same volume as the constant rectangular beam. 

This may be accomplished by assuming a constant width 
for all designs, integrating the power law (4.1) and setting 
the resulting value equal to that of the classical rectang- 
ular beam. Using symmetry and integrating the power law 
l 

from x=0 to x= ^ 



1/2 
f h 
0 



(1 - ^) n 
o^ a ’ 



dx = V 



(4.29) 



determines the volume of the variable design. For the rec- 
tangular beam the volume for one-half the structure is 
h r bd 

— — . Thus, equating the two determines 

h Q = (n+1 ) h r (4.30) 



It is obvious that the center height must be expanded for 

the variable design in order that the volumes are compatible. 

When h> is restricted to a fixed value the critical 
o 

loads may be calculated as described in Section IV-A, i.e. 



39 



using the intersecting value of (4.21a) and (4.21b) 

P£ 2 



Recalling that K = 



8 






j the critical load is defined as 



8(2-n)x N /c°C°' 



crit 



(4.31) 



where x represents the intersecting value of the Bessel 
equation curves. The amended load is now found by sub- 
stituting for C° C° and h . 

12 o 



crit 



8 ( 2-n) xh 

,2 Vb 6 EG/36 . 



(4.32) 



Thus j 



8(2-n)x(n+l)h 

P* = — 3 (4.33) 

where 

c 3 ■ IF (4-34) 

Cj, = EG. (4.35) 

The amended load is referred to as P* since its value is 
greater than its unexpanded counterpart. It is left in 
the above form in order that it may be readily compared to 
the solution for the rectangular beam. 

A computer solution is included in the Computer Program 
Section and provides for determining the intersecting value 
of (4.21a) and (4.21b). It additionally calculates the 
critical and amended loads for each value of considered n. 
These results are summarized in Table I below. 
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TABLE I 



n 


V 


v+1 


Pcrit 


P* 


P*/Pu 


0 


0.2500 


1.2500 


16 . 936=Pu 


16.936 


1.0 


.1 


0.2368 


1.2368 


15.618 


17.179 


1.014 


.2 


0.2222 


1.2222 


14.288 


17.145 


1.012 


• 3 


0.2059 


1.2059 


12.944 


16.827 


0.994 


.4 


0.1875 


1.1875 


11.580 


16.212 


0.957 


.5 


0.1667 


1.1667 


IO.I89 


15.283 


0.902 


.6 


0.1429 


1.1429 


8.758 


14.013 


0.827 


.7 


0.1154 


1.1154 


7.263 


12.347 


0.729 


.8 


0.0833 


1.0833 


5.656 


10.181 


0.601 


.9 


0.0455 


1.0455 


3.794 


7.210 


0.426 


0.98 


0.0010 


1.0010 


0.507 


1.013 


0.059 


1.0 


0 


— 


0 


0 


0 



The best design exists between n=0.1 and 0.2. Additional 
calculations indicated that n=0.l4 offers the largest buckling 
load (P*=17.199) for the variable height formulation. This 
design is shown in Figure 8. Calculation of values for 
P*/Pu indicates that there is little change for the first six 
values of n. A graph at P# versus n is shown in Figure 9- 
Figures 27-38 in Appendix E illustrate the various designs 
studied in the interval of n=0 to n=l. 
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Y-SCRLE,-.-2.0GE : .+00 UNITS INCH. 

ROBERT L. BURNS L LDR. USN, THESIS 
URRIRBLE HEIGHT BERM DESIGN — CASE N=0.14 

Figure 8. Simply Supported Beam Design - 

Case h = h (1- 2x/L) 0 ' 1^ P* = 17.199. 




N 

Figure 9- Amended Load P* Versus Values of N, (case 1). 

C. EXTENSION OF THE PROBLEM 

The formulation for design criteria was previously based 
on the assumption 

h = h Q (l - 2x/£) n . (4.1) 

The case will now be considered for the design based on the 
condition 

h = h Q (l + 2x/£) n . (4.36) 

This condition is physically shown in Figure 10. By inspec- 
tion of the governing equations of the forementioned solution, 
it is evident that the sign change makes no difference in 



4 3 




Figure 10. Variable Height Beam Based on h = h Q (l+2x/£) n . 



the solution of the problem. The identical critical loads 
are again found. However, each design is restricted to the 
same volume of material. This requires that the midspan 
height be decreased while the end heights of the beam are 
increased. This modification has a significant effect on 
the amended loads . 

Integrating (4.36) and using the same technique described 
in Section IV-B, the center height reduces to 



n+1 



( 2 ) 



n+1 



- 1 



(4.37) 



while the ends expand to 



jn+l 

h e ~ 2 (n+1 ) h r 



(4.38) 



A computer solution for this problem is included at the end 
of the thesis. The results are summarized in Table II. 

The structure is less stable for any comparable condi- 
tion listed in the previous section. Figures 39-48, 
Appendix E include designs analyzed for this problem. A 
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TABLE II 



n 


X 


Pcrit 


P* 


% Decrease 


0 


1.05851 


16.936 


16.936 


— 


.1 


1.02747 


15.618 


15.023 


11.30 


.2 


0.99220 


14.288 


13.215 


21.97 


.3 


0.95178 


12.944 


11.508 


32.05 


.4 


0.90470 


11.580 


9.891 


41.60 


.5 


0.84907 


10.189 


8.359 


50.64 


.6 


0.78200 


8.758 


6.898 


59.27 


.7 


0.69840 


7.263 


5.490 


67.58 


.8 


0.58920 


5.656 


4.102 


75.78 


.9 


0.43119 


3.794 


2.639 


84.41 



plot of the amended load P* versus n is shown in Figure 11. 
Note that all designs reflect a decrease in buckling load 
as compared to Figure 9 and Table I. Also, there is no 
range of designs for which the value P* remains relatively 
stationary . 



D. LATERAL BUCKLING OF A CANTILEVER 

A similar analysis [7] for a cantilever beam shows that 
the buckling- loads are determined by 



J -v ( 2^> ■ 0 



where K is defined as 



K = 



pr 
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(4.39) 



(4.40) 




N 

Figure 11. Amended Load P* Versus Values of N (case 2). 

A computer solution was used to calculate the first zero 
of a negative order Eessel function of the first kind, and 
used this result to obtain the critical and amended loads 
for the structure. These results are summarized below in 
Table III. The last value is obtained by evaluation of the 
Bessel function of the second kind and use of tables [ 1 8 ] . 
The values of P* versus n are shown in Figure 12. Note in 
particular that all values are greater than the ultimate 
critical load, Pu, for the rectangular case. 

The best design is determined for the case of n=0.75 5 
and is noted as P*=M.9l8. This shape is shown in Figure 13* 
The range of shapes considered in the formulation are shown 
in Appendix E. 
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TABLE III 



n 


V 


Pcrit 


P# 


% Increase 


0 


-.2500 


4 . 013=Pu 


4.013 


— 


.1 


-.2368 


3.858 


4.239 


5.63 


.2 


-.2222 


3.694 


4.433 


10.47 


.3 


-.2059 


3.534 


4.595 


14.50 


.4 


-.1875 


3.374 


4 . 724 


17.72 


.5 


-.1667 


3.214 . 


4 . 821 


20.13 


. 6 


-.1429 


3.053 


4.885 


21.73 


.7 


-.1154 


2.892 


4.917 


22.53 


.8 


-.0833 


2.731 


4.915 


22.48 


.9 


-.0455 


2.568 


4.880 


21.60 


1.0 


0 


2.405- 


4 . 810 


19.86 



E. COMPARATIVE ANALYSIS OF SOLUTIONS 

A comparative analysis between the solutions determined 
for the cantilever and simply supported beams provides some 
interesting insight into the actual physics of the problem. 
The results obtained from the cantilever solution appear to 
be a near optimum shape for the structure. At first sight 
analysis of the simply supported beam does not appear to 
yield the optimum shape, although a better design than 
that of classical solution is obtained. However, a close 
examination of the physics of the problem leads to a predic- 
tion for the possible optimum shape, although it was not 
analytically determined. 
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Figure 12. Amended Load P* Versus Values of N 
for Cantilever Eeam. 



The value n=0.l4 provided a best buckling load of 
P* = 17.199 for the simply supported beam. This corresponds 
to a gain of only 1.55? over that of a constant rectangular 
cross section. For the cantilever, it was found that 
P* = 4.918 (n=0 . 75 ) which is a 22.5 % gain over that of the 
rectangular section. Physically these results can be 
explained when the lateral bending and twisting moments of 

the two problems are examined. 

Recalling equations (3-4) and (3-6) for the rectangular 

beam yields the following principle axes moments 
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( 4 . 42 ) 




Evaluation at the boundary condition, x=0 : u^=u and u'=0 

shows that the twisting moment M^=0 and bending moment 
assumes its maximum value. At the ends of the structure, 
x = 2'} <f> = 0 and u=0. This implies that bending moment M^=0 
and twisting moment is a maximum. 

The combination of these two moments appear to restrict 
the beam from assuming an optimum shape which differs 
significantly from the uniform case. In fact, if the 
structures depicted in Appendix E are examined, it may be 
noted that as the beam approaches a pyramidal shape, the 
less stable the structure becomes. Physically, this seems 
to imply that the vanishing end area offers less and less 
resistance against the twisting moment, and the structure 
will readily buckle. 

It is apparent and expected that the existing moment 
conditions govern the optimum shape of the structure. 
Physical intuition would indicate that the maximum amount 
of material be located where the maximum moment existed. 

If and were equal and a linear variation could be 
assumed, the resulting optimum shape would actually be 
that of a rectangular beam. For the problem considered, 
appears to have a domineering influence on the optimum 
design. If this was not the case, then the solutions 
determined in Section IV-C would have indicated a "best" 
buckling load greater than those for the present solution. 
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The actual moment values for M and M r are not known 
and cannot be determined unless the displacement and twist 
functions u and <j> are known exactly. Additionally, these 
functions are coupled to the variable cross sectional area 
as indicated by equations (3-3 a,b,c). Analyzing equations 
(4.4l) and (4.42) and Figure 4 leads to the following inter- 
pretation. The moment representation of appears to be 
parabolic while that of M (based on a parabolic <f>) appears 
to be cubic and exhibits an inflection point. These patterns 

are shown in Figure 14. The maximum combined moment must 

£ 

then occur between x=0 and x= ^ • 



M^ = lateral bending moment 




Figure 14. Expected Moment Configuration and 
Approximate Optimum Design. 



On the basis of moment evaluation alone, the optimum 
design would apparently be of a form similar to the shape 
shown in Figure 14. Since is expected to be greater 
than Mj. the design would lump the largest amount of material 

l 

somewhere between x=0 and x= ^ but closest to x=0. 
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Available time has prevented analysis of this expected 
"optimum" shape. A Rayleigh Ritz formulation for this 
design was attempted by assuming the displacement and twist 
functions and using a Lagrange interpolating polynomial to 
approximate the design. The results indicated a buckling 
load slightly less than that of the uniform rectangular 
case . 

Because of the low gain in efficiency (1.55$) from the 
best design obtained from the variable height solution, it 
appears that the true optimum shape will not reflect a 
large performance gain over that of the rectangular case. 

In fact, the determined value from the variable height 
solution is expected to be very close to the actual optimum 
design . 

P. VALUES OF N GREATER THAN ONE 

The analyses of the problems discussed have been restric- 
ted to the range for n=l. If n is greater than one the 
actual curvature of the beam design becomes concave. For 
the simply supported beam this physically implies that the 
structure must become less stable than the limiting pyramidal 
shape. This is evident since even less area is distributed 
toward the beam ends. But, it was found earlier that the 
pyramidal shape is actually unstable for any load. By this 
reasoning the problem must be restricted to the range of 
values used to have physical meaning. 

Mathematically, if the value of n is allowed to be 
greater than one, then equation (*J.15) changes to the form 
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<p = y 



( 4 . 45 ) 



1-n 

2 [a 1 J_ v (p) + a 2 J v (p) ] 

for which must be zero. The resulting roots introduce 
a singularity condition and the solutions are misrepresen- 
tative of the problem. Consequently, it is deemed necessary 
for the value of n to be restricted to the interval 
0 < n <_ 1 to provide meaningful results. 
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V. CONCLUSIONS 



For the simply supported beam the assumed power law 
restricts the number of designs which may be analyzed. 
However, a good indication of the problem behavior can be 
determined from the results. The number of designs which 
remain relatively stable do not differ radically from the 
uniform case. As the designs converge to a pyramidal 
shape they rapidly become less stable. 

It appears that the combination of the two moments , 
and govern the design. As the beam approaches a 
pyramidal shape, less area becomes available to resist the 
twisting moment and the beam readily buckles. Thus the 
interaction of the two moments restrict the beam from 
attaining a shape which differs significantly from the 
rectangular case. 

The effect of the moment values on beam design is 
better demonstrated from the cantilever results. At the 
free end of the beam there is no moment and the design calls 
for very little material to be located at that position. 

The interaction of the two moments seem to produce a maximum 
slightly away from the fixed end - otherwise, the design 
for the best buckling load would have indicated the pyra- 
midal shape. The cantilever solution, thus, helps to 
support the validity of the simply supported beam results. 



5 4 



The true optimum shape for the simply supported beam 
has not been determined by the variable height formulation. 
However, the expected optimum shape was approximated from 
analyzing the physics of the problem. This shape is not 
expected to offer any significant increase over the design 
obtained from the variable height formulation. 
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APPENDIX A 



A STRUCTURAL OPTIMIZATION PROBLEM 

An example of a structural optimization problem which 
illustrates the usefulness of the variational method of 
Section II is that of obtaining the optimum shape of a uni- 
formly loaded beam. This problem is a "strength" problem 
in contrast to the buckling problem which exhibits no 
strength. For the problem considered, bending is in the 
vertical plane of symmetry (Fig. 15). The height of the 




Figure 15- Variable Width Uniformly Loaded Beam. 

beam is taken as constant, and a uniform load is applied. 

1 . Statement of the Problem. 

The design b(x) is to be determined for a variable 
width which provides a maximum strength capacity for the 
structure. Unless some specific constraints are imposed on 
the problem, results may well be meaningless. For instance, 
if a minimum volume structure was considered with no other 
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